First-order framework and generalized global defect solutions 
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This work deals with defect structures in models described by scalar fields. The investigations 
focus on generalized models, with the kinetic term modified to allow for a diversity of possibilities. 
We develop a new framework, in which we search for first-order differential equations which solve 
the equations of motion. The main issue concerns the introduction of a new function, which works 
like the superpotential usually considered in the standard situation. We investigate the problem in 
the general case, with an arbitrary number of fields, and we present several explicit examples in the 
case of a single real scalar field. 

PACS numbers: 11.27.+d 



I. Introduction. The first recent observations col- 
lected in [l( have led us with the intriguing fact that the 
Universe is presently undertaking accelerated expansion. 
These informations directly contributed to establish some 
important advances in Cosmology, one of them being the 
presence of dark energy, which is introduced to respond 
for more then 2/3 of the contents of the Universe - see, 
e.g., Refs. [2J for some recent reviews of the subject. 

The presence of dark energy has opened some distinct 
routes of investigations. The general idea starts with 
the recognization that the standard cosmological model 
is the ACDM model, which describes gravity with the 
standard scalar curvature R, with the inclusion of the 
cosmological constant and non relativistic dust-like mat- 
ter. The dark energy modification of this scenario may 
be considered with the inclusion of several distinct pos- 
sibilities, the two most active directions of investigations 
being the quintessence way, in which one considers the 
possibility of changing the standard scenario with the in- 
clusion of scalar fields, and the f(R) way, in which one 
modifies gravity, changing the contribution of the scalar 
curvature R to f(R), with f(R) being a non trivial func- 
tion of the scalar curvature. 

In the present work we shall follow the quintessence 
way, that is, we shall deal with standard gravity but 
changing the scalar field dynamics. The use of non canon- 
ical models has appeared both in early time inflation and 
in dark energyscenarios to describe the current phase of 
the Universe [3|, [j, Q. But here, however, we shall make 
an important reduction which concerns the absence of ge- 
ometrical degrees of freedom. That is, we shall deal with 
flat space-time. The idea of freezing out the geometri- 
cal degrees of freedom is strategical: we shall first deal 
with general scalar field dynamics, and then, in another 
work, we shall pursue the yet more general scenario, ac- 
tivating the geometrical degrees of freedom [6[ . The final 
goal is to investigate the general model, in which we mix 
quintessence and f(R), but here we first deal with the 
simpler case, where the scalar fields are allowed to inter- 
act and self-interact in the flat space-time. 

Such investigations are directly related to the recent 
investigations 0, @, QJ, LL0, LLU, LL2| where one discusses 
general properties of topological defects such as domain 



wall, strings and monopoles [13J. These investigations 
have dealt with formal aspects of unidimensional topo- 
logical solutions 7:] and specific study on kinks and the 
corresponding linear stability [8j, including thepresence 
of compactons for a specific class of models Q. More 
recently the case of local vortex-like defects have been 
investigated in [lOj] and similar modifications with higher- 
order scalar kinetic terms in the braneworld scenario with 
a single extra dimension [ll|, [l2| • 

Although the presence of new non linear effects which 
come from the generalized dynamics complicates the 
search for analytical solutions, we take motivation from 
the previous investigations [Tj| to consider the possibil- 
ity of finding first-order differential equations that solve 
the corresponding equations of motion. This possibility 
is implemented with the introduction of a general func- 
tion, W, which in principle depends on all the scalar fields 
and nicely leads to the presence of first-order equations. 
As usual, the first-order equations simplify the investi- 
gation as we show below, where we split the subject in 
two parts, one dealing mainly with the formal aspects of 
defects and the other, which is dedicated to the study of 
specific examples. 

II. Generalities. We start with a model given by 
a non standard action, described by the set of n scalar 
fields {<pi,4>2, ■■■,4>n} in the two-dimensional space-time 



S = / d 2 x C(<f>i,Xjk) 



(1) 



where i,j,k — 1,2, ...,n. We are considering the metric 
(H — ), and the scalar fields are all dimensionless, with 
dimensionless space and time coordinates, and coupling 
constants. Here the quantities Xij are defined by 



X-; 



1 



• d^id^ 



>Vj 



(2) 



This is the action which controls the generalized mod- 
els. It allows obtaining the energy-momentum tensor in 
the form 



Tpv = Lx^d^idyipj - g^C 



(3) 



and we are using the notation: C a = dC/da, that is, 
£x tj — dC/dXij and C^ = dC/dcpi, etc. We shall be 



searching for static solutions, and for static fields we can 
write the energy and stress densities, p(x) — T 00 and 
t(x) = Tn, as respectively 



p(x) 



■C 



r(x) = Cx^W + C 



(4a) 
(4b) 



The equations of motion for the n real scalar fields 
fa — fa{x, t) have the non standard form 



d fJt (Lx ii d^ j )=C 4>i 

We expand these equations to get 



Gfd a 



■2Xji£, Xii <i>i - £4 



where 



Gg = C Xti ri afl + £x a x im 8 a <, 



IV <p m 



(5) 



(6) 



(7) 



For static solutions fa — fa(x), the equations of motion 
reduce to 



C 



A", 



or better, 



C^Xij + 2CxuX im Xlm)<t > j = ^XjiCxij^i — £4 



(8) 



(9) 



These equations can be integrated to give L — 
2Cxi Xij = C, where C is an integration constant which 
can be identified as the stress density (jib]) . Stability of 
the static solution requires the vanishing of C, that is, 
stability requires the stressless condition |8J. Therefore 
we write 



£ — 2Cx i:l Xij 







(10) 



Notice that this equation depends of the scalar fields and 
their first derivatives. Therefore, it is a first-order equa- 
tion, but it is a constraint equation, which we name the 
stressless constraint. 

We now deal with the search of first-order differential 
equations. We consider the energy density, which for the 
stressless solutions can be written as 



p = -£ = C XiJ < 



(11) 



In this case, if we introduce the new function W — 
W(fa,fa>,-;fai) such that 

C Xt] fa) = W^ (12) 

we obtain the energy density in the form 



,-*„«-£ 



(13) 



in a way such that the energy can be written as the vari- 
ation of W, that is 

E = AW = w{M<x>),<h(po),...,M<x>)) 

- W{M-oo),M-oo),...,M-oo))(U) 



In order to circumvent problem with unstable solutions, 
here we assume that W is a C 1 function. 

We substitute (fl"2")) in the equations of motion @ to 
get 



W A 



(15) 



which is also a set of first-order differential equations 
which solve the equations of motion. Since the set of 
equations (fl"2"|) is a first-order set of equation, we name 
the procedure the first-order framework for the scalar 
field models under investigation. In this sense, the above 
investigation falls into the program shown in [lj|, in 
which we dealt with a diversity of models, searching for 
first-order differential equations which solve the corre- 
sponding equations of motion. The first-order equations 
are also important to investigate linear stability, because 
the appearance of W eases the general investigation and 
helps to factorize the Schrodingcr-like Hamiltonian into 
the two first-order differential operators S and S' , which 
we explicitly construct below. 

The above calculations generalize the standard sce- 
nario, which appears when the scalar fields evolve under 
usual dynamics. We note that the energy does not de- 
pend on the explicit form of the solution, but it is given 
in terms of W, calculated at the asymptotic values of 
the fields. The non vanishing of the energy indicates the 
presence of non trivial field configuration having non triv- 
ial topology which ensures existence and stability of the 
solution. 

The stability of the static solution can be inferred by 
the presence of the topological current, defined in terms 
of the function W introduced above: 



f T = £ ^d v W 



(16) 



This definition is inspired in [15j , and it is constructed to 
make the topological charge directly related to the energy 
of the static solution. We see that for static solution 
jj, = dW/dx, so that Qt is directly related to AW. 

We now focus on linear stability. We consider 
fa(x,t) = fa (x) +r]i (x, t), supposing that rji(x, t) are small 
fluctuations around the static solution. In this case we 
get, going up to first-order in the fluctuations, 



with 



Xi 



1 



*J ' ij 



1 



Xij = -d^fad^ + -d^jd 11 ^ 



(17) 
(18) 



The procedure gives 

£<j>i —> £<pi + £4>i4>jVj + £-<t>iX jk Xjk (19a) 

C Xlj -> C Xi] + £x t jfaVk +Cx il x mj Xi m (19b) 

We use these expressions into the equations of motion (JSJ) 
to get 

d^Cxv&TH + C Xt] x ]m d^faX lm ) 
= [C^ 3 - d^Cx^d^fan)] Vj (20) 



In the case of static solutions one gets 

C x tJ Vj ~ {(£x it +C Xil x im X lm )r]' j )' 

= {£<n<t>j + i^x im 4> 3 (f> m Y)v] ( 21 ) 

We suppose that the fluctuations are given by 

r)i(x,t) = r)i(x) cos(uit) (22) 

in order to get 

- [(jC Xii +Cx a x im X lm )^]' 

= [C-fabj + (Cxi^j'Pm)' +u 2 C Xi] ] Vj (23) 

This equation has the general form 

(-aijV'jY = Mi (24) 

We can modify the above equation into the Schrodinger- 
like equation 

/ ,72 \ 

(25) 



5ij — -j + Uij ) uj — w u 



where the potential U is now a matrix which depends on 
the matrix S and R, introduced as folows: in (j!M|) we 
change r]i(x) into Ui(z), such that 



Tjl OijUj 

dz 

ax = — 

R 



(26a) 
(26b) 



In this case, the Schrodinger-like equation requires that 
t dSjk d(a,ijR) 



2aijR , 

dz dz 



S jk = (27a) 

/ "u,^x m] = kj (27b) 



R-^S^aTlC 



The general study is awkward. The issue here is that it 
is better to investigate linear stability in a case by case 
manner. We will return to this below, where we examine 
the given examples specifically. 

The above investigations are done on general grounds. 
To see how it generalizes the standard situation, let us 
suppose that the Lagrange density has the very specific 
form 



C = M lJ X i:j - V 



(28) 



where M y = M y (0i) and in principle M may also de- 
pend on the scalar fields. In this case we have that, using 
the first-order equation (Ql 



(29) 



b< = M^W^ 



The stressless constraint (TIUl) leads to 



V = ^W^M-^W^ 



(30) 



which shows that the function W is a generalization of 
the superpotential which appears in the standard sce- 
nario. To be explicit, let us consider the standard case, 
in which M is constant, M y = 5 % K Here we have 



C = ^d^id^j - V(<P) 



Also, from (fl2|) we can write 
and from (1 101) we have to have 



V(<P) 



ISijWfrWi 



(31) 



(32) 



(33) 



This is the standard scenario and this the reason why we 
use W to represent the new function that we had to in- 
troduce in (|12|L As we know, in the standard case we can 
include fermions to get to the supersymmetric extension 
of the model, and there W is named the superpotential. 
This poses a nice issue, which concerns the inclusion of 
fermions in the case of generalized dynamics, to get to 
the supersymmetric model - this is an interesting issue 
which is presently under consideration. 

The very interesting feature of the present generaliza- 
tion is that we get to first-order differential equations 
which solves the equations of motion, even though we do 
not specify the potential. This is a strong result, which 
shows that we can get to the first-order framework even 
though we do not know explicitly how the scalar fields 
interact formally. 

III. Applications. Let us now move toward applica- 
tions. For simplicity, we consider models described by a 
single real scalar field, with the equation of motion being 
written in terms of W. The general model is described by 



C = C{X, 



(34) 



This case was already investigated in [3, H, Q under spe- 
cific conditions, but here we will introduce other improve- 
ments. The equation of motion and the stressless condi- 
tions are given by 



and 



[£*</>']' 



C - 2£ X X = 



The first-order equation has the form 
Cx<t>' = Ws 



(35) 



(36) 



(37) 



This is the general case. The first-order equation sup- 
ports kinklike solutions, and can have the topological 
charge shown before. The study of linear stability follows 
the general investigation given in the former section. We 
solve the equations (f2"6")l and (|2"7| to get 



dx 



dz 
T 



and r\ 



JCx~A 



(38) 



where A 2 — (2£xxX + Cx)/ &x- This allows writing 
the Schrodinger-like equation 



— u zz + U(z)u = UJ U 



(39) 



where 



U(z) 



(VA£x~) z 
VAC 



x 



1 

Tx~ 



£00 + -j ( £<px-^ 



(40) 



See also Ref. [8j for other details. 

We illustrate the general behavior with the model 



C = X - aX 2 - V{4>) 



or equivalently 



C = X + aX\X\~V((f>) 



(41) 



(42) 



for a > real and positive parameter. In this case, the 
equation of motion is 



(1 - 3a(j)'' 



= V rh 



(43) 



and the first-order and stressless equations are given by, 
respectively 



(t,' + a^ = Wa, 



3q 



Vtf) 



(44a) 
(44b) 



The first equation is an algebraic equation of third degree 
in <f>' . The only real solution is 



<t>' = G a (W+) 
where the function G a (W^) is such that 

9a(v) 2 



G a (v) 



Go 



9a{v) 



(45) 



(46) 



and 



g a (v) = (54a 2 v + 6V3(16 a 3 + 27a 4 v 2 ) 1/2 ) 1/3 (47) 

We now use the stressless constraint to write the poten- 
tial in the form 



v{<t>) = \cl m) 



3q 



Gtm) 



(48) 



The advantage of the introduction of the function W 
is that we know a priori the energy of the topological so- 
lutions, since their asymptotic values are found from the 
relation G l f,(W ( j > ) = 0. In this case, V(4>) is a-dependent. 
The inclusion of the quartic term changes some of the 
features of the solutions, as the width and energy. We 
illustrate this with the function 



W 



1 



We use this W to plot in Fig. 
and energy density 



1 the potential, solution 



P (x) = G 2 a (W^)+aG 4 a (W tP ) 



(50) 



which are depicted for some specific values of the param- 
eter a, with a = being the thicker line. The first-order 
equation <\>' = G a (l — a 2 ) is solved numerically. In Fig. 1, 
we see that both the potential, and the width and energy 
density of the solution change with a, showing how the 
modification introduced by a change the specific features 
of the defect solutions. 




FIG. 1: Plots of the potential (upper panel), kinklike solution 
(middle panel) and energy density (lower panel) for the model 
with W — 4>— 4> 3 /3, with several values of a. The case a = is 
shown with the solid line, and the cases with a — 0.5, 1.0, 2.0, 
appear with dashed, dot-dashed and dotted line, respectively. 



The configurations have asymptotic values at —1 for 
x — > — oo, and +1 for x — > oo. Therefore, as the 
energy only depends on the asymptotic values, E = 
W(4>(+oo)) — W{4>{— oo)), the solution have the same 
energy of the standard cf) A case, which is E = 4/3. In 
Ref. |8j, we have studied models in which the potential 
had no explicit dependence on a, as well as the solution. 
v^JJ Here, in the present models we see that it is the energy 



that does not change with a. In this sense, the present 
investigation adds to the former scenario of [|| another 
interesting possibility, in which the energy of the topo- 
logical solution is fixed, although its specific form may 
change together with the potential. 

If we consider a very small we can go further and write, 
up to first order in a 



G a (v) = v — av 3 



(51) 



which is good approximation for v limited. We note that 
Gq(v) = v and G a (0) = 0, and in this case the equations 
(1431) lead to 



4>' = Wt- aWl 
and so the potential becomes 

We then integrate (|52")) to get to 

4>{x) = 4> Q {x) - aW^Wifo) 



(52) 



(53) 



(54) 



where 4>o( x ) is the solution for a — > 0. The corresponding 
energy density is 

p(x) = p (x) (1 - a(Wl + 2WW Mo )) (55) 

where p (x) = Wj o . 

The explicit form W = <fi — (j> 3 /3 leads to the standard 
results 



4>o(x) = tanh(a;) 
Po{x) = sech (x) 



(56a) 
(56b) 



which are valid for a equal to zero; for a very small, we 
get the new analytical results up to first order in a 

0(x)=tanh(a;)-^tanh(x)sech 2 (x)(2 + sech 2 (x))(57a) 
p(x) = sech 4 (x) (l - ^ (7 sech 4 (x) + 4 sech 2 (x) - 8)V57b) 

There is no contribution to the energy from the 
a— dependent terms in p in (I57bj) . and this is in accor- 
dance with the fact that the energy is E = AW = 4/3 
and only depends on W and the asymptotic values of the 
field, which do not depend on a, as we see from (|57al) . 
We use 



dx = (l + aW$)dz 

r, = (1-aWhu 



(58a) 
(58b) 



in order to find the potential of the Schr6dingcr-like equa- 
tion 



U(z) = ^% + Wl 



(59) 



and for W given by (|49]) we get 

U(z) = 4- 6 sech 2 (z) 
2 



1 



a (sech 4 (z) - 2 sech 2 (z) + 1) j (60) 



and the zero mode is given by, explicitly 



uq(z) = -sech 2 (z) 



1 - -a (5 sech 4 (z) - secli (z) - 1) ) (61) 



We can consider another model, given by 

C = X\X\-V(4>) (62) 

or yet, we make it more general and consider the case 



£ = 



-X\X\ 



- v{4>) 



(63) 



where n = 1,2, ... is non vanishing, positive integer. The 
case n = 1 leads to the standard model, and for n = 2 
we get to dHll), which is the model recently investigated 
in HP. 

The equation of motion for n arbitrary is given by 



dV 



(64) 



{2n-l)^ n - l U" = — 
d(j) 

We use the first-order equation to write 

0' - Wf=* (65) 



and the stressless condition leads to the potential 

In- 1 



V{4>) 



2a 



Wi 



and we also have 

2n 

p = (j) = w^ 
We also have A 2 = In — 1 and the quantity 

V^A-(2n-l) 1 /4 2 ^|0'|«-i 



(66) 



(67) 



(68) 



Since we are searching for topological solutions, the 
scalar field is supposed to be monotonic, so we take ip' > 
0. Also, we change from dx — ► dz and r\ — > u in order to 
get 



dx 



1 



-dz 



s/2n - 1 

77 = ^n-l)- 1 ^^^') 1 - 



(69a) 
(69b) 



and so the Schrodinger-like equation can be written in 
the form 



(70) 



Hu(z) = ( -— + U(z) ) u{z) = w^u(z) 



where the potential is 



fV(d>) 



r)(n — 9\ 4(n-l) 2.1-3 

U(z) = -^TT^ 2 "" W^+nW^W^ (71) 

We factorize this equation with the first-order operators 
5 and St such that H = S^S, with 



5 



2(r»-l) 

= w k=i W, 



dz V2n- 1 
In this case the normalized zero mode is given by 



(72) 



(73) 



where E is the energy of the corresponding solution. 

We now consider the case n = 2 explicitly. Here the 
equation of motion becomes 



U 



12 ill 



dV 



The first-order equation is 



6' = W 1 ' 3 



(74) 



(75) 



and the stressless condition leads to the potential 

V(<P) = TWt /3 (76) 



We also have A 2 — 3, and the fluctuations about 
the kinklike solution with </>' non negative leads to the 
Schrodinger-like equation 

Hu{z) = f-|j + ZW-^W^A u(z) = w 2 u(z) (77) 

It can be factorized with the first-order operators S and 
S+, with 



d 2 V3 2 /3 TJ/ 



(78) 



This result shows that there is no bound state with neg- 
ative eigenvalue, and the normalized zero mode is given 
by 



u (z) 



W. 



2/3 



(79) 



We can further illustrate the above case with the spe- 
cific functions 



W 



5 r V 



(80a) 



W = ^v/l - 02 PL _ I^\ + 3 arcsin((/)) (gob) 



which give the following potentials 

3 
4 
3 
4 



V{4>) 



(i-0 2 ) 4 
(1-0 2 ) 2 



(81a) 

(81b) 




FIG. 2: Plots of the potentials of the models given by (|81a[) 
(solid line) and ()81b|l (dashed line). 



respectively, which we plot in Fig. [2]. 

In the first case the model will lead us with kinklikc 
solution. In the second case the field behaves differently, 
and the model will lead us with compactons, in a way 
similar to the case considered in [9j. 

Let us now study the model described by (|80a[) . The 
first-order equation is cf>' — 1 — <fi 2 and we get the kinklike 
solution <p(x) — tanh(a;), with energy E — AW — 32/35. 
For the fluctuations, the potential in the Schrodinger-like 
equation gets the form 



U{z) = 48-60sech 2 (v / 3z) 



(82) 



This is the modified Poschl- Teller potential [16| and the 
normalized zero mode is now 



uq(z) 



1 



70v / 3sech 4 (v / 3~ z ) 



(83) 



In Fig. 3 we plot both U(z) and the zero mode uq{z). 
This potential supports the zero mode and other three 
bound states, with energies E r = 3r(8 — r), for r — 
0,1,2,3. 




FIG. 3: Plots of the potential (solid line) and zero mode 
(dashed line) given by (|82[) and (|83|l . respectively. 



The other model is described by (|80b[) . The first-order 
equation is cf>' = (1 — cj) 2 ) 1 ' 2 , and may support compactons 
[l7J, which are growing in importance in a diversity of 
scenarios, in particular from the point of view of pattern 
formation, since patterns usually appear in nature with 
finite extent. In the present work, the explicit solution 
of the first-order equation is 



x < -§ 



0, for 
c) = < sin(x), for 

0, for x > f 



o — X — 2 



(84) 



The corresponding energy is E = 3tt/8. In this case the 
potential and the normalized zero mode are given by 

oo, for z < -jj^j 
U{z)={ -12 + 6sec 2 (V3z), for - -S- < z < -S- (85) 



°°> for Z> 171 



and 



uq(z)-- 



, 0, for z < --=7= 

93 2V3 

cos 2 (V3z), for - -5= < z < -5= 



7r%/3 



°. fOT ^>^ 



2V3 — — 2V3 



The potential in (|85p is the Poschl- Teller potential [lq . 
which is plotted in Fig. 4, together with the above zero 
mode. 




In this case, U{z) has an hybrid asymptotic behavior, 
having Poschl- Teller and modified Poschl- Teller tails, as 
we show in Fig. 6. The energy of the hybrid configuration 
is E = 1024v / 2/H55, and the zero mode has the form 



u (z) 



^J= tanh 2 h§z\ sech 4 (fz),forz>0 

0, for z < 

(91) 
which is also plotted in Fig. 6. In this case, the reflection 
coefficient is unit, and this would modify the standard 
behavior concerning the scattering of particles by the wall 
itself, changing the frictional force which would act on the 
wall, introducing new features in the wall evolution in a 
(3,1) cosmological scenario [13j . but this is out of the 
scope of this work. 




FIG. 5: Plot of the potential given by (JS7J). 



FIG. 4: Plots of the potential (solid line) and zero mode 
(dashed line) given by (|85p and (|86p . respectively. 

The interesting feature of the potential in ([55)) is that 
it only supports bound states, and for r = 0, 1, 2, ..., the 
corresponding eigenvalues are given by E r — 12r(r + 2). 
Its behavior is similar to the potentials found in Q for 
models of the form C = V((f>)F(X), which appears in the 
case of tachyons and their generalizations. 

The case of compactons may lead to an hybrid config- 
uration, named half-compacton, as noted in [9] for the 
model with the potential 



' U(z) 



v^) = -(i-^ni+^r 



(87) 



which is plotted in Fig. 5. The solution has the form 

ri-2tanh 2 (^),forz>0 
I 1, for x < 



Here we get 

2 



W 



1155 



(533 + 755</>-455</> 2 + lO5(/> 3 )(l-0) 5/2 (89) 



The potential for the Schrodinger-like equation is such 
that 



8-28 sech 2 (^A +21 scch 4 (-#z) 



U(z) = { tanh 2 (^2 

oo, for z < 



for z > 



(90) 




FIG. 6: Plots of the potential (solid line) and zero mode 
(dashed line) given by (|90p and (|9ip . respectively. 

IV. Ending comments. In this work we have shown 
how to extend the first-order formalism to generalized 
models, described by scalar fields with non standard dy- 
namics. The main issue here is the introduction of W and 
the related first-order equations of motion (fT2")) . These 
equations together with the stressless condition (TT5|) sim- 
plify the investigations of linearly stable defect solutions, 
which in specific conditions can be obtained analytically. 

We have investigated several models, in general de- 
scribed by a set of real scalar fields. In the case of a single 
field, we have dealt with the Lagrange density of the form 
F(X) — V((f>), and we have considered the cases where 
F(X) is governed by X + aX\X\ or X\X\ n -\ We have 
also considered several explicit forms for W , in particular 
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the case with W — <f> — 3 /3, which reproduces the stan- 
dard </> 4 model if the kinetic term is also of the standard 
form F(X) = X. The general result is that topological 
defects also appear in models with generalized dynam- 
ics, and there it is also possible to formulate a first-order 
framework in which one deals with first-order differen- 
tial equations, which solve the corresponding equations 
of motion and ease the study concerning the search for 
analytical results. 

In the present work, we have focused mainly on ex- 
plicit models governed by a single real scalar field. Ev- 
idently, we can also deal with two or more fields, and 
we will return to the subject in a future work, where we 
investigate more sophisticated models, governed by two 
or more real or complex scalar fields. In particular, we 
will deal with generalized models constructed from com- 



plex fields ipi , in which the complex fields have the form 
fi = 4>i + iXi- I n this case, the function F(X) should 
deal with X — d t Jp i d'- i (pi 1 with the sum over all the com- 
plex fields implied. This form is specific but important, 
since it allows generalizing the dynamics one usually finds 
in the standard Wess-Zumino model, which has interest- 
ing features which will certainly help us to make explicit 
integrations [LSI ]. 

The presence of topological defects in models described 
by scalar fields is well known in high energy physics [1 31 ] , 
and the present work contributes to show how they ap- 
pear in models where the dynamics is generalized to in- 
clude higher order contributions on X. The results are 
of current interest to high energy physics, since nowa- 
days there are situations in which modifications of the 
standard dynamics are welcome. 
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